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Abstract 

Campisi, Zhan, Talkner, and Hanggi state, in promoting a new logarithmic computational ther- 
mostat [ arXiv 1203.5968 and 1204.4412 ] , that (thermostated) Nose-Hoover mechanics is not 
Hamiltonian. First I point out that Dettmann clearly showed the Hamiltonian nature of Nose- 
Hoover mechanics. The trajectories { q(t) } generated by Dettmann's Hamiltonian are identical to 
those generated by Nose-Hoover mechanics. I also observe that when the (Hamiltonian) Campisi 
thermostat is applied to "nonequilibrium" heat transfer problems some very interesting, and some- 
what paradoxical, phase portraits result. See too Marc Melendez' nice arXiv 1205.3478 as well as 
our joint work arXiv 1206.0188. 
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In 1984 Shuichi Nose discovered a deterministic, time-reversible, logarithmic thermostat 1 . 
This computational thermostat imposes a time- averaged kinetic temperature kT = m( v 2 ) 
through a thermostated Hamiltonian. His Hamiltonian includes a "time-scaling variable", 
s, along with its conjugate momentum p s . Consider the simplest interesting example. For a 
single harmonic oscillator (with mass m, force constant k, relaxation time r and Boltzmann's 
constant k all set equal to unity) Nose's thermostated Hamiltonian is : 

2-Hnosc = [ (p/s) 2 + q 2 +p 2 s + Tln(s 2 ) ] . 

The equations of motion which follow from Nose's Hamiltonian , 

q = (p/ s 2 ) ; p=-q; s = p s ; p s = (p 2 /s 3 ) - (T/s) , 

are somewhat "stiff" because s can become arbitrarily small. "Scaling the time" in these 
equations of motion, by a factor of s, gives a new, and better behaved, set : 

q = (p/s) ; p = -sq ; s = sp s ; p s = (p/s) 2 - T . 

Then, to convert to the simpler "Nose-Hoover" form: introduce v = (p/s) and ( = p s : 

q = v ; q = v = -q-^q=- q -( V] ( = [ q 2 - T } = [ v 2 - T } . 

The time scaling used to obtain the Nose-Hoover equations suggests (wrongly, it turns 
out) that they are "non-Hamiltonian" . Though well-behaved, these Nose- Hoover motion 
equations are not ergodic. For the harmonic oscillator they have a wide variety of periodic, 
nearly-periodic, and chaotic, solutions 2 . 

On the other hand, in July 1996, Carl Dettmann discovered that a different, but closely 
related, Hamiltonian , 

2%Dettmann = S"H Nos 6 = (p 2 / ' s) + s[ q 2 + p\ + Thl(s 2 ) ] = , 

gives directly the Nose-Hoover motion equations but without any time scaling 5 . The equa- 
tions of motion from Dettmann's Hamiltonian are as follows : 

q = (p/s) ■ p = - sq ; s = sp s ; p s = (l/2)(p/s) 2 - [ (l/2)q 2 + (l/2)p 2 s + Tins ] - T . 

Because Dettmann's Hamiltonian is identically equal to zero, the combination in square 
brackets is equal to —(1/2) (p/s) 2 : 

[ (l/2)g 2 + (l/2)(p 2 + Tins ] = ~(l/2)(p/s) 2 . 
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Introducing v = (p/s) and ( = p s again produces the Nose-Hoover equations, but this time 
without the need for any time scaling : 

q = v ; q = v = -q-^q=-q-( v - ( = [ q 2 - T ] = [ v 2 - T ] . 

A year later Dettmann and Gary Morriss published yet another Hamiltonian forrn^ : 

2H DM = e- Q [ p 2 + P 2 ] + e +Q [ q 2 + 2TQ } = . 

This Dettmann-Morriss Hamiltonian, though slightly different, leads, in a very similar way, 
to the same Nose-Hoover motion equations. This approach was soon rediscovered by Bond, 
Leimkuhler, and Laird 6 , who refer to Dettmann's original Hamiltonian as the Nose-Poincare 
Hamiltonian. All these discoveries make the Hamiltonian nature of the Nose-Hoover equa- 
tions quite clear. 

Campisi et alii 1 cited Klages' book 8 as the source of their "nonHamiltonian" characteriza- 
tion of Nose- Hoover mechanics. In that book Klages warns the reader: "one should carefully 
distinguish between traditional Hamiltonian formulations of classical mechanics and the gen- 
eralized Hamiltonian formalism outlined [ in the book ]." It seems to me that one can take 
this warning too seriously. In several helpful and stimulating private communications Camp- 
isi, Hanggi, and Klages emphasized the arbitrary value of Dettmann's Hamiltonian, zero, as 
well as the failure of the (q, v, () Nose-Hoover equations to follow Liouville's Theorem. 

Admittedly, Dettmann's zeroing of the Hamiltonian could be accomplished in several 
ways. One way is to add a term proportional to s , chosen in conjunction with the initial 
conditions in order to make the Hamiltonian vanish. This is essentially the "Poincare" trans- 
formation noticed by Bond, Leimkuhler, and Laird. Continuing with the simple harmonic 
oscillator example, with T additionally chosen equal to unity : 

{ q, p, s,p s } = { 1, 0, 1, 1 } [ initial conditions ] , 

the Hamiltonian becomes : 

^Another = (s/2) [ {p/s) 2 + q 2 + p\ } + s ln(s) - s . 

The corresponding equations of motion that follow from it (and satisfy Liouville's theorem) 
are : 

{ Q = (P/s) ; P = -sq ; s = sp s ; p s = (l/2)[ (p / s) 2 - q 2 - p 2 s ] - Ins } . 
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The solution of these equations is identical to that from the Nose-Hoover equations, 

{q = v ; i) = -q-£v\ C = v 2 - 1 } , 

where v = (p/s), ( = p s , and the initial values are { q, v, ( } = { 1, 0, 1 } . 

There is no obvious way to introduce a second temperature into the Nose, Dettmann, or 
Dettmann-Morriss Hamiltonian. Nevertheless multi-temperature problems can be treated 
easily by generalizing the Nose-Hoover equations of motion to control sets of velocities 
through a set of friction coefficients^. By sandwiching Newtonian degrees of freedom between 
two sets of boundary particles (a "hot" set and a "cold" set) it is easy to simulate steady heat 
flow. It is well-established that such dissipative systems give multifractal strange attractors 
in the full system+thermostats phase spaced 

Campisi et alii claim that their own Hamiltonian , 

Hcampisi = H usual + (kT/2) ln(5 2 + S 2 ) + (P 2 /2M) , 

plus an unspecified (and crucial) weak coupling between the system variables { q } and the 
thermostat variable 5* is an improvement^. If their thermostat is more easily matched in 
laboratory experiments then it is indeed a step forward. But, if a multi-temperature Campisi 
Hamiltonian, including £[ (kTi/2) ln(<5 2 + Sf) ], could impose more than a single temper- 
ature on selected degrees of freedom, either Liouville's theorem or the fractal phase-space 
structures that arise away from equilibrium would be casualties. Typically, Hamiltonians 
don't give fractals. 
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6x6 {q,p} Phase-Plane Plots 



Figure 1: Phase-plane portraits with e = 0.01 . The specified boundary temperatures are 0.5 and 
1.5 . Actual kinetic temperatures are { 0.46,0.70,0.74, 1.40 }, left to right, starting at the bottom. 

I. ADDENDUM OF 30 APRIL 2012 

How can a strictly Hamiltonian thermostat satisfy Liouville's Theorem (incompressible 
flow in phase space) while simultaneously giving the multifractal structures (on all scales) 
associated with dissipative systems? These two demands are contradictory, and so merit 
investigation. After trying several inconclusive modifications of the simple thermostated 
harmonic oscillator problem I introduced more complexity, considering a short (two-particle) 
"0 4 " chain, with one particle coupled to a "hot" Campisi Thermostat and the other to 
a "cold" one. The Figures show phase portraits for both the chain and the thermostat 
particles, for the two choices {e} = { 0.01,0.10 }, using the Hamiltonian : 

2U = j2lTMS* + e) + Pn + 4(S 1 -x 1 ) 2 + (S 2 -x 2 ) 2 ] + (x 4 1 + xl)/2 + (x 1 -x 2 ) 2 + pj+pi. 
i=i 

The particle phase-plane { q,p } trajectories are shown in Figures 1 and 2. The two imposed 
temperatures, nominally 1.5 and 0.5, and the weak-coupling parameter e = 0.01 and 0.10 , 
provide typical, but quite surprising, results. 

The nature of these solutions is at the least "odd", from the standpoint of the Second 
Law of Thermodynamics. They provide a large temperature gradient, but with nowhere 
for the heat to flow! There are no energy sources or sinks in these problems. Evidently 
these Hamiltonian systems cannot be dissipative, despite their temperature gradients, and 
have instead perfectly well-behaved solutions satisfying Liouville's incompressible Theorem, 
(df/dt) = . 

This odd behavior deserves a more thorough investigation than I can undertake here. 
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Weak Coupling 
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8x8 {q,p} Phase-Plane Plots 



Figure 2: Phase-plane portraits with e = 0.10 . The specified boundary temperatures are 0.5 and 
1.5 . Actual kinetic temperatures are { 0.71, 0.53, 0.60, 0.97 }, left to right, starting at the bottom. 

But the broad outlines seem clear. The nature of the Hamiltonian thermostat coupling is 
crucial. If the coupling is weak then the thermostats oscillate independently of the "system" 
(here a two-particle chain). If the coupling is strong then the thermostats no longer impose 
the desired kinetic temperatures on the hot and cold parts of the system. These problems 
deserve, and will no doubt receive, detailed investigation. 
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II. ADDENDUM OF 23 MAY 2012 



There are two other examples of non-dissipative Hamiltonian thermostats in "Hamilto- 
nian Dynamics of Thermostated Systems: Two-Temperature Heat-Conducting 4 Chains" , 
published in the Journal of Chemical Physics 126, 164113 (2007). The "Hoover-Leete Hamil- 
tonian thermostat" keeps the kinetic energy constant through a straightforward application 
of Lagrangian constraints. The "Landau-Lifshitz Hamiltonian thermostat" keeps the config- 
urational temperature constant. Both these latter thermostats are unable to stimulate heat 
flow. The log-thermostat of Campisi et alii [ Reference 7, which will soon appear in Physical 
Review Letters ] is certainly the simplest example of this anomalous behavior. 

I would particularly like to thank Michele Campisi, Peter Hanggi, Rainer Klages, and 
Marc Melendez for a series of stimulating and educational emails on this subject. See also 
Marc's very recent and perceptive [ ar^iv 1205.3478 ] contribution of 15 May 2012: "On the 
Logarithmic Oscillator as a Thermostat" . 

III. ADDENDUM OF 16 JUNE 2012 

More details of the log-thermostats' shortcomings will appear in Marc Melendez' and my 
joint work, arXiv 1206.0188 . Another example of Hamiltonian bases for "nonHamiltonian" 
thermostated equations of motion appears in Wm. G. Hoover, B. Moran, C. G. Hoover, and 
W. J. Evans' "Irreversibility in the Galton Board via Conservative Classical and Quantum 
Hamiltonian and Gaussian Dynamics", Physics Letters A 133, 114-120 (1988) . 
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